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Calibration,  Verification,  and  Application  of  a 
Two-Dimensional  Flow  Model 


D.  Michael  Gee*  M.ASCE 


Abstract 


Information  was  required  concerning  velocity  distributions  in  the  outlet 
channel  immediately  downstream  from  the  Harry  S.  Truman  Dam  generating 
facility  (see  Fig  1) .  This  information  was  to  be  used  to  ascertain 
hydraulic  forces  and  flow  directions  for  structural  design  of  a  fish  net 
or  other  type  of  fish  barrier  to  be  placed  across  the  outlet  channel. 

The  velocity  distribution  in  the  outlet  channel  is  complex,  governed  by 
operation  of  various  combinations  of  1  to  6  pump-turbine  units  in  either 
generation  or  pumpback  mode.  The  flow  pattern  can  be  further  complicated 
by  spillway  releases  with  or  without  simultaneous  operation  of  the  power¬ 
house.  To  provide  timely  design  information  for  the  planned  fish  facili¬ 
ties,  a  mathematical  model  was  used  to  predict  the  flow  fields. 

Selection  of  a  Mathematical  Model 


A  model  for  simulation  of  two-dimensional  free  surface  flows  in  the  hori¬ 
zontal  plane  was  deemed  appropriate  for  this  study.  The  problem  was 
particularly  well  suited  for  use  of  the  two-dimensional  finite  element 
hydrodynamics  model,  RMA-2?p  which  had  been  used  previously  by  The  Hydro- 
logic  Engineering  Center  (HEC)  on  several  project  applications l1 * 

Data  Summary 

The  data  required  to  perform  this  study  may  be  divided  into  three  cate¬ 
gories:  (1)  "Run"  data,  i.e.  that  information  required  to  execute  a 
simulation,  (2)  "Calibration"  data,  i.e.  prototype  measurements  which 
are  used  to  adjust  various  model  coefficients  to  bring  the  model's  per¬ 
formance  into  conformance  with  that  of  the  prototype,  and  (3)  "Verifi¬ 
cation"  data;  additional  prototype  measurements  used  to  evaluate  model 
performance . 

Run  data  were  derived  from  construction  drawings  of  the  outlet  channel 
and  related  physical  features.  Calibration  and  verification  data  con¬ 
sisted  of  several  sets  of  detailed  velocity  measurements.  Both  magnitude 
and  direction  of  flow  velocity  were  measured  at  several  points  in  the 
vertical  at  several  locations  across  the  channel  (see  Fig  1)  for  three 
discharges  as  shown  in  Table  1.  Values  used  in  this  study  were  averages 
of  the  point  vertical  data  at  each  location. 


*  Research  Hydraulic  Engineer,  Hydrologic  Engineering  Center 
US  Army  Corps  of  Engineers,  Davis,  California,  USA 

Presented  at  the  ASCE  Conference  on  Frontiers  in  Hydraulic  Engineering  held 
at  the  Massachusetts  Institute  of  Technology,  Cambridge,  Mass,  on  9-12  Aug¬ 
ust  1983. 


TABLE  1.  CALIBRATION  CONDITIONS 


DATE 

DISCHARGE 

TAILWATER  ELEV 

UNITS 

MEASURED 

(cfs) 

(ft.  NCVD) 

5  May  82 

27,000 

656.0 

4,5,6 

26  May  82 

24,900 

658.8 

4,5,6 

27  May  82 

15,280 

658.1 

4  &  5 

Note:  1  cfs 

=  0.028  m3/s;  1  ft 

=  0.305  m. 

Development  of 

the  Finite  Element 

Network/Problem  Schematization 

The  finite  element  network  was  developed  to  provide  detail  in  the  area 
where  a  fish  net  was  being  considered  and  in  areas  of  anticipated  strong 
velocity  gradients.  Sufficient  network  detail  was  provided  at  the  power¬ 
house  face  to  resolve  flows  emanating  from  individual  units.  The  area 
schematized  extended  from  the  powerhouse  face  downstream  approximately 
2700  feet. 

The  spillway  stilling  basin  was  also  included  in  the  area  modeled  as 
circulations  in  that  area  are  important  to  the  flow  field  in  general  and 
to  allow  for  simulation  of  spillway  flows  if  desired.  In  general, 
curved-sided  elements  were  used  along  the  flow  boundaries  to  allow  tan¬ 
gential  flow  along  the  banks.  The  network  used  for  calibration  and  pro¬ 
duction  runs  is  shown  on  Fig.  1.  Note  the  distinct  rows  of  elements 
along  the  possible  net  alignments;  this  provides  a  mechanism  for  simula¬ 
ting  the  effect  of  the  fish  net  on  the  flow  field  by  increasing  the  bed 
roughness  within  one  of  these  rows  of  elements  to  account  for  head  loss 
across  the  net.  Boundary  conditions  used  were;  inflow  (generate)  or  out¬ 
flow  (pumpback)  rates  at  the  appropriate  units,  slip  conditions  for  other 
flow  boundaries,  and  a  prescribed  water  surface  elevation  at  the  down¬ 
stream  boundary.  All  simulations  were  of  steady  state  conditions. 

Calibration  Process/Results 

The  data  obtained  on  26  May  82  were  used  to  calibrate  the  model  and  the 
other  two  data  sets  used  for  verification.  The  measured  velocities 
showed  a  flow  reversal  on  the  left  side  of  the  channel  with  velocities 
there  of  up  to  1  ft/sec  (0.305  m/sec)  directed  towards  the  powerhouse. 
This  circulation  is  apparently  driven  by  flows  from  the  powerhouse  per¬ 
sisting  as  a  jet  along  the  right  bank.  Initial  model  runs  indicated  a 
tendency  for  the  jet  to  stay  on  the  right  side  of  the  channel,  however, 
it  diffused  and  mixed  across  the  channel  much  more  completely  than  was 
observed  in  the  prototype.  The  simulated  velocities  were  in  the  down¬ 
stream  direction  across  nearly  the  entire  cross  section  at  the  location 
of  the  measurement.  Use  of  smaller  turbulent  exchange  coefficients  was 
indicated  to  decrease  mixing  in  the  model.  However,  when  coefficients 
less  than  about  15/ft2/sec  (1.4m2/sec)  were  used  the  Newton-Raphson 
method  used  to  solve  the  nonlinear  system  of  equations  did  not  converge. 

Boundary  roughness  is  described  by  the  Chezy  equation  in  RMA-2.  A  Chezy 
coefficient  of  120  ft'^/sec  (66  nri/sec)  was  used;  approximately  equivalent 
to  a  Manning's  n  of  0.020.  Bottom  friction  does  not  play  a  dominant  role 
in  this  problem,  as  evidenced  by  the  small  head  loss  in  the  reach  of 


interest  (this  was  confirmed  by  preliminary  sensitivity  runs).  Improve¬ 
ment  of  the  calculated  flow  distribution  by  varying  the  bottom  roughness 
spatially  vas  not  attempted  because  there  was  no  physical  justification 
for  doing  so. 

The  initial  approach  to  obtaining  a  convergent  solution  with  lower  tur¬ 
bulent  exchange  coefficients  was  to  increase  network  detail.  (The 
initial  network  contained  228  elements  and  509  nodes;  the  final  313  ele¬ 
ments  and  686  nodes.)  It  was  thought  that  this  would  improve  the  model's 
performance  for  two  reasons:  (1)  increased  network  detail  would  allow 
resolution  of  smaller  scale  flow  features,  and  (2)  smaller  exchange  co¬ 
efficients  are  usually,  used  with  smaller  elements.  Model  performance 
improved  somewhat,  but  was  still  unsatisfactory.  Attempts  were  also  made 
to  vary  the  turbulent  exchange  coefficients  in  the  transverse  direction 
from  those  in  the  longitudinal  direction.  These  efforts  again  did  not 
significantly  improve  the  results. 

Convergence  of  the  Newton-Raphson  algorithm  requires  that  the  first 
"guess"  at  the  solution  be  fairly  close  to  the  final  solution;  even  for 
relatively  well-behaved  functions.  The  first  guess  used  in  RMA-2  is 
zero  velocity  everywhere  and  a  horizontal  water  surface.  Provision  also 
exists  to  use  as  the  first  guess  the  solution  from  a  previously  run  simu¬ 
lation.  HEC  had  used  this  feature  in  a  previous  study  of  a  situation 
involving  a  large  change  in  water  surface  elevation  within  the  study 
domain.  Successive  problems  were  solved,  each  having  larger  water  sur¬ 
face  gradients  and  providing  the  first  guess  at  the  subsequent  problem 
solution  until  the  desired  problem  was  solved.  It  was  decided  to  try 
this  technique  by  using  the  convergent  solution  obtained  with  large  ex¬ 
change  coefficients  as  the  beginning  guess  for  a  solution  with  a  smaller 
exchange  coefficients.  This  technique  proved  successful  in  allowing  use 
of  exchange  coefficients  of  7.7  ft2/sec  (0.7m2/sec)  which  gave  a  solu¬ 
tion  that  closely  matched  the  observed  data  (see  Fig  2). 

The  other  critical  aspect  of  successfully  modeling  the  H.S.  Truman  flow 
field  was  selection  of  the  depth  at  the  powerhouse.  The  prototype  dis¬ 
charge  enters  (or  exits)  over  about  the  lower  2/3  of  the  tailrace  depth 
and  is  not  uniformly  distributed  over  this  portion.  As  RMA-2  is  a  verti¬ 
cally  averaged  model,  capturing  this  vertical  detail  was  not  possible. 
Consequently,  proper  entrance  velocities  were  obtained  by  adjusting  the 
bottom  elevation  at,  and  immediately  downstream  from,  the  powerhouse  face. 

The  model  was  then  tested  against  the  verification  data  sets  (see  for 
example  Fig.  3)  and  produced  results  equally  as  good  as  the  calibration 
run . 

Production  Runs 


Once  the  model  was  calibrated  and  verified,  production  runs  were  made. 

All  coefficients  and  physical  characteristics  of  the  model  were  unchanged 
from  the  calibration  runs;  only  inflow  locations  and  quantities  and  down¬ 
stream  tailwater  elevations  were  changed  to  model  the  various  operational 
configurations.  To  date,  fourteen  generation  and  pumpback  configurations 
have  been  analyzed  with  the  model.  The  computed  velocity  fields  have 
provided  valuable  information  for  the  selection  and  design  of  a  viable 
fish  facility. 
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Fig.  4  Data  Flow  and  Program  Linkage 
Computational  Aspects 

Performance  of  this  type  of  study  requires  manipulation  of  several  data 
files,  application  of  at  least  three  different  computer  programs,  and 
graphics  capability  for  analysis  of  input  data  and  simulation  results. 

The  data  flow  and  correspondence  between  the  computer  programs  are  shown 
on  Fig.  4.  All  simulations  were  performed  on  HF.C '  s  Harris  500  mini¬ 
computer.  The  RMA-2  simulations  required  3-4  minutes  of  c.p.u.  time 
per  iteration,  or  a  total  time  of  about  20  minutes  per  simulation. 
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C.  P.  Davis  and  A.  J.  Eredrlch, 

Aug  71,  20  pp . 

TP-32  Alternative  Approaches  to  Water  Resource 
System  Simulation,  L.  R.  Beard, 

A.  Weiss,  and  T.  A.  Austin,  May  72, 

12  pp. 

TP-33  System  Simulation  for  Integrated  Use 

of  Hydroelectric  and  Thermal  Power 
Generation,  A.  J.  Eredrlch  and 
L .  R.  Beard,  Oct  72,  22  pp. 

TP-34  Optimizing  Flood  Control  Allocation  for  a 
Multipurpose  Reservoir,  f.  K.  Duren 
and  L.  R.  Beard,  Aug  72,  15  pp. 

TP-35  Computer  Models  for  Ralnfa 1 1 -Runoff  and  River  ADA-104  909 

Hydraulic  Analysis,  D.  W.  Davis, 

Mar  73,  46  pp. 

TP-36  Evaluation  of  Drought  Effects  at  Lake  At  1 1 1  an , 

A.  D.  Feldman,  Sep  72,  15  pp. 

This  publication  Is  not  available  to 
countries  outside  of  the  U.S. 

TP-37  Downstream  Effects  of  the  Levee  Overtopping  at  ADA-104  908 

Wilkes-Barre,  PA,  During  Tropical  Storm 
Agnes,  A.  D.  Feldman,  Apr  73,  22  pp. 

TP-38  Water  Quality  Evaluation  of  Aquatic  Systems,  ADA-104  907 

R.  G.  Willey,  Apr  73,  22  pp. 

TP-39  A  Method  for  Analyzing  Effects  of  Dam  ADA-104  906 

Failures  In  Design  Studies, 

W.  A.  Thorns,  Aug  72,  29  pp. 

TP-40  Storm  Drainage  and  Urban  Region  Flood  Control  ADA-104  920 

Planning,  0.  W.  Davis,  Oct  74,  40  pp. 


ADA-104  912 


ADA-104  911 


ADA-104  910 


V 


TECHNICAL  PAPERS  (TP)( Continued ) 


HE  C 

NUMBER 

TP-41 

TP-42 

1 P  -43 

1 P  -44 

TP-46 

TP-46 

TP-47 

TP-48 

1  P-49 

TP-50 


TITLE 


HEC  NT1S 

PRICE  NUMBER 


>2.00  Each 


HEC-5C,  A  Simulation  Model  for  System 
formulation  and  Evaluation, 

B.  S.  Elchert,  Mar  74.  28  pp. 

Optimal  Sizing  of  Urban  Flood  Control 

Systems,  0.  W.  Oavls,  Mar  74,  18  pp. 

Hydrologic  and  Economic  Simulation  of  flood 
Control  Aspects  of  Water  Resources 
Systems,  B.  S.  Elchert,  Aug  75,  10  pp. 

Sizing  Flood  Control  Reservoir  Systems  by 
Systems  Analysis,  B.  S.  Elchert  and 
0.  W.  Davis,  Mar  76,  34  pp . 

Techniques  for  Real-Time  Operation  of  flood 
Control  Reservoirs  In  the  Merrimack 
River  Basin,  B.  S.  Elchert,  J.  C.  Peters 
and  A.  f.  Pabst.  Nov  75,  45  pp. 

Spatial  Data  Analysis  of  Nonstructural 

Measures,  R.  P.  Webb  and  M.  W.  Burnham, 
Aug  76,  21  pp. 

Comprehensive  flood  Plain  Studies  Using 
Spatial  Data  Management  Techniques, 

D.  W.  Davis,  Oct  76,  20  pp. 

Direct  Runoff  Hydrograph  Parameters  Versus 
Urbanization,  D.  L.  Gundlach, 

Sep  76,  7  pp. 

Experience  of  HEC  In  Disseminating 

Information  on  Hydrological  Models, 

B.  S.  Elchert,  Jun  77,  9  pp. 

(superseded  by  TP#56). 

Effects  of  0am  Removal:  An  Approach  to 
Sedimentation,  0.  T.  Williams, 

Oct  77,  36  pp. 


ADA-106  242 

ADA-106  243 

ADA-106  244 

ADA-106  245 

ADA-106  246 

ADA-106  247 

ADA-106  248 

ADA-106  249 


ADA-106  250 


TECHNICAL  PAPERS  ( TP H Continued) 


HEC  HEC  NTIS 

NUMBER  TITLE _ _ PRICE  NUMBER 


j2 . 00  Each 

TP-51  Design  of  Flood  Control  Improvements  by  ADA-106  364 

Systems  Analysis:  A  Case  Study, 

H.  0.  Reese,  A.  V.  Robbins, 

J.  R.  Jordan,  and  H.  V.  Ooyal, 

Oct  71 ,  23  pp. 

TP-52  Potential  Use  of  Digital  Computer  Ground  ADA-106  251 

Water  Models,  0.  L.  Gundlach, 

Apr  78,  38  pp . 

TP-53  Development  of  Generalized  Free  Surface  ADA-106  252 

Flow  Models  Using  Finite  Element 
Techniques,  0.  M.  Gee  and 
R.  C.  MacArthur.  Jul  78,  21  pp. 

TP-54  Adjustment  of  Peak  Discharge  Rates  for  ADA-106  253 

Urbanization,  0.  L.  Gundlach, 

Sep  78,  7  pp. 

TP-55  The  Development  and  Servicing  of  Spatial  ADA-106  254 

Data  Management  Techniques  In  the 
Corps  of  Engineers,  R.  P.  Webb  and 
D.  W.  Davis,  Jul  78,  26  pp. 

TP-56  Experiences  of  the  Hydrologic  Engineering  ADA-106  255 

Center  In  Maintaining  Widely  Used 
Hydrologic  and  Water  Resource 
Computer  Models,  B.  S.  El  chert, 

Nov  78,  16  pp. 

TP-57  Flood  Damage  Assessments  Using  Spatial  ADA-106  256 

Data  Management  Techniques,  D.  W.  Davis 
and  R.  P.  Webb,  May  78,  27  pp. 

TP-58  A  Model  for  Evaluating  Runoff -Quality  In  ADA-106  257 

Metropolitan  Master  Planning, 

L.  A.  Roesner,  H.  M.  Nlchandros, 

R.  P.  Shublnskl,  A.  D.  Feldman, 

J.  W.  Abbott,  and  A.  0.  Frledland, 

Apr  72,  81  pp. 


i 

f 


HE  C 

NUMBER 

TP-59 

TP-60 

TP-feT 

TP-62 

TP-63 

TP-64 

TP-65 

TP-66 

TP-67 


TECHNICAL  PAPERS  ( TP ) ( Continued ) 


TITLE 


HEC  NTIS 

PRICE  NUMBER 


>2.00  Each 


Testing  of  Several  Runoff  Models  on  an 
Urban  Watershed,  J.  Abbott, 

Oct  78,  53  pp. 

Operational  Simulation  of  a  Reservoir 
System  with  Pumped  Storage, 

G.  f.  McMahon,  V.  R,  Bonner  and 
B.  S.  Elchert,  Feb  79,  32  pp. 

Technical  Factors  In  Small  Hydropower 
Planning,  D.  W.  Davis,  Feb  79, 

35  pp. 

Flood  Hydrograph  and  Peak  Flow  Frequency 
Analysis.  A.  D.  Feldman.  Mar  79  21  pp. 

HEC  Contribution  to  Reservoir  System 
Operation,  B.  S.  Elchert  and 
V.  R.  Bonner,  Aug  79,  28  pp. 

Determining  Peak-Olscharge  Frequencies  In 
an  Urbanizing  Watershed:  A  Case  Study, 

S.  F.  Daly  and  J.  C.  Peters,  Jul  79,  15  pp. 

Feasibility  Analysis  In  Small  Hydropower 
Planning,  D.  W.  Davis  and  B.  W.  Smith, 

Aug  79.  20  pp. 

Reservoir  Storage  Determination  by  Computer 
Simulation  of  Flood  Control  and 
Conservation  Systems,  B.  S.  Elchert, 

Oct  79,  10  pp. 

Hydrologic  Land  Use  CTassIf Icatlon  Using 
LANOSAT,  R.  J.  Cermak,  A.  D.  Feldman 
and  R.  P.  Webb,  Oct  79,  26  pp. 


ADA-106  258 

ADA-106  259 

ADA-109  757 

ADA-109  758 

ADA-109  759 

ADA-109  760 

ADA-109  761 

ADA-109  762 

ADA-109  763 


TP-68  Interactive  Nonstructural  Flood-Control 
Plannng,  D.  T.  Ford,  Tun  80,  12  pp. 


ADA-109  764 


HEC 

NUMBER 

TP-69 

TP-70 

TP-71 

TP-72 

TP-73 

TP-74 

TP-75 

TP-76 

TP-77 

TP-70 


TECHNICAL  PAPERS  ( TP  HCont  Inued ) 


TITLE 


HEC  NTIS 

PRICE  NUMBER 


>2.00  Each 


Critical  Hater  Surface  by  Minimum  Specific  ADA-951  599 

Energy  Using  the  Parabolic  Method, 

B.  S.  Elchert,  1969,  15  pp . 

Corps  of  Engineers  Experience  with  ADA-109  765 

Automatic  Calibration  of  a 
Precipitation-Runoff  Model,  D.  T.  Ford, 

E.  C.  Morris,  and  A.  D.  Feldman. 

May  80,  12  pp. 

Determination  of  Land  Use  from  Satellite  ADA-109  766 

Imagery  for  Input  to  Hydrologic  Models. 

R.  P.  Webb,  R.  Cermak,  and  A.  D.  Feldman, 

Apr  80.  18  pp. 

Application  of  the  Finite  Element  Method  to  ADA-109  767 

Vertically  Stratified  Hydrodynamic  Flow 
and  Hater  Quality,  R.  C.  MacArthur  and 
U.  R.  Norton,  May  80,  12  pp. 

Flood  Mitigation  Planning  Using  HEC-SAM,  ADA-109  756 

D.  W.  Davis,  Jun  80,  17  pp. 

Hydrographs  by  Single  Linear  Reservoir  ADA-109  768 

Model,  J.  T.  Pederson,  J.  C.  Peters, 
and  0.  J.  Helweg,  May  80,  17  pp. 


HEC  Activities  In  Reservoir  Analysis, 

V.  R.  Bonner,  Jun  80,  10  pp. 

Institutional  Support  of  Water  Resource 
Models,  J.  C.  Peters,  May  80,  23  pp. 

Investigation  of  Soil  Conservation  Service 
Urban  Hydrology  Techniques, 

D.  G.  Altman,  W.  H.  Espey,  Jr.  and 
A.  0.  Feldman,  May  80,  14  pp. 

Potential  for  Increasing  the  Output  of 
Existing  Hydroelectric  Plants, 

D.  W.  Davis  and  J.  J.  Buckley, 

Jun  81 ,  20  pp . 


ADA-109  769 

ADA-109  770 

ADA-109  771 


ADA-109  772 


6 


HEC 

NUMBER 

TP-79 

TP -80 

TP-81 

TP-8? 

IP-83 

TP-84 

TP-85 

TP-86 

TP-87 

TP -88 

TP -89 


TECHNICAL  PAPERS  ( TP H Continued) 


TITLE 


HEC  NTIS 

PRICE  NUMBER 


*2.00  Each 


Potential  Energy  and  Capacity  Gains  from 
flood  Control  Storage  Reallocation 
at  Existing  U.  S.  Hydropower 
Reservoirs,  B.  S.  Elchert  and 

V.  R.  Bonner,  Jun  81.  18  pp. 

Use  of  Non-Sequentlal  Techniques  In  the 

Analysis  of  Power  Potential  at  Storage 
Projects,  G.  M.  franc,  Jun  81,  18  pp. 

Data  Management  Systems  for  Water  Resources 
Planning.  D.  W.  Davis,  Aug  81,  12  pp. 

The  New  HEC-1  Mood  Hydrograph  Package,  A.  D. 
feldman,  P.  B.  Ely  and  D.  M.  Goldman, 

May  81 ,  28  pp . 

River  and  Reservoir  Systems  Water  Quality 
Modeling  Capability,  R.  G.  Willey, 

Apr  82,  15  pp. 

Generalized  Real-Time  flood  Control  System 
Model,  B.  S.  Elchert  and  A.  F.  Pabst, 

Apr  82,  18  pp. 

Operation  Policy  Analysis:  Sam  Rayburn 
Reservoir,  0.  T.  ford,  R.  Garland 
and  C.  Sullivan,  Oct  81,  16  pp. 

Training  the  Practitioner:  The  Hydrologic 
Engineering  Center  Program, 

W.  K.  Johnson,  Oct  81,  20  pp. 

Documentation  Needs  for  Water  Resources 

Models,  W.  K.  Johnson,  Aug  82,  16  pp. 

Reservoir  System  Regulation  for  Water 
Quality  Control,  R.G.  Willey, 

Mar  83,  18  pp. 

A  Software  System  to  Aid  In  Making  Real-Time 
Water  Control  Decisions,  A.  f.  Pabst 
and  J.  C.  Peters,  Sep  83,  17  pp. 


ADA-109  787 

ADA-109  788 

ADA-114  650 

ADA-114  360 

ADA- 1 1 4  19 2 

ADA-114  359 

ADA-123  526 

ADA-123  568 

ADA-123  558 

ADA- 130  829 


? 


TECHNICAL  PAPERS  (TPHContlnued) 


TITLE 


[ 


i 


Calibration,  Verification 
of  a  Two-Dimensional  Flow 
D.  Michael  Gee.  Sep  83,  6 


and  Application 
Model . 

PP- 


